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ABSTRACT 

A new method for measuring gravitational lensing with high redshift type la supernovae 
f— ^ ■ is investigated. The method utihzes correlations between foreground galaxies and supernova 

C^ ! brightnesses to substantially reduce possible systematic errors and increase the signal to noise 

ratio. It is shown that this lensing signal can be related to the mass, size and substructure of 
galaxy halos and galaxy clusters if dark matter consists of microscopic particles. This technique 
may be particularly useful for measuring the size of dark matter halos, a measurement to 
which the lensing of galaxies is not well suited, and for measuring the level of substructure in 
galaxy halos, a problematic prediction of the cold dark matter model. The contributions to the 
signal from galaxy halos and galaxy clusters are modeled and contributions to the noise from 
04 ' fluctuations in the galaxy number counts, galaxy redshift error, dispersion in SN luminosities 

J^ , and sample variance are estimated. The intrinsic distribution of supernova luminosities and 

Qs I its redshift evolution are removed as major sources of uncertainty. The method is found to be 

complimentary to galaxy-galaxy lensing. The required observations of ^ 100 supernovae have 
already been proposed for the purposes of cosmological parameter estimation. 
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H ' 1- Introduction 

^ ' Outside of the visible extent of galaxies not a great deal is known about the distribution of dark 

matter (DM) on scales smaller than galaxy clusters. Galaxy-galaxy lensing (Fischer et al. 1999; Hudson 
k> ; et al. 1998; Brainerd, Blandford, & Small 1996; Griffiths et al. 1996; Tyson et al. 1984) and studies of 

$H ' the satellite galaxies orbiting the Milky Way and other galaxies (Wilkinson & Evans 1999; Kochanek 1996; 

Zaritsky & White 1994) have put some constraints on the mass of galaxy halos, but their size scale and the 
degree to which they are smooth mass distributions or collections of subclumps are not well determined. 
Even the concept of galaxy halos, with a single galaxy within each of them and well defined sizes, may 
not be correct. Nor is it clear how the observable properties of galaxies relate to the DM around them. 
It would be preferably if these properties could be derived directly from cosmological initial conditions 
and the properties of DM, but this has proven to be an extremely difficult problem. Simulations of galaxy 
formation give incomplete answers which are hard to compare with observations. For example pure cold 
dark matter (CDM) simulations contain a large amount of small scale structure which is not mirrored in the 
observed distribution of light (Moore et al. 1999a). The process of galaxy formation is highly dependent 
on cooling and feedback from stars which must be modeled ad hoc. As a result complete galaxies have 
never been created in a hydrodynamic simulation. In addition, changing the properties of DM so that it 
is self-interacting or a mixture of hot and cold DM, for example, can have a strong influences on small 
scale structure. The theory of galaxy formation, and DM, would clearly benefit from further observational 
constraints on these scales. Gravitational lensing is presently the most promising method available for 
probing the distribution of dark matter and its relationship with light. 



The discovery of a tight correlation between the hght curve shape and the luminosities of type la 
supernovae (SNe) has resulted in their great utility as cosmological distance indicators (Hamuy et al. 1996; 
Riess, Press, & Kirshner 1996) and as a possible probe of gravitational magnification. The measured 
standard deviation of corrected peak luminosities in local SNe is now ^ 0.12 mag. Two groups have used 
these SNe to constrain the distance-redshift relation and through this cosmological parameters (Perlmutter 
et al. 1997; Riess et al. 1998). The successes of these studies have inspired plans for more aggressive, larger 
searches for high redshift SNe in the near future. The volume and quality of data is likely to increase 
dramatically. These data will have implications well beyond constraints on f2m and Qa - the average mass 
density and vacuum energy density. Other uses include measuring the universe's equation of state and star 
formation history. Here I discuss the gravitational lensing of high redshift SNe and its implications for dark 
matter, its distribution and its correlation with light. 

Metcalf & Silk (1999) have already discussed using the lensing of SNe to determining the composition of 
DM. It was found that DM made predominantly of macroscopic compact objects with masses ^ lO^'^ Mq 
(ex. MACHOs, primordial black holes) can be differentiated from DM made of microscopic particles (ex. 
WIMPS, axions) through studying the magnification probability distribution of ~ 100 SNe at z ~ 1. In 
the case of macroscopic DM, the most likely brightness of a supernova (SN) is well below the mean value, 
and near the empty beam solution or Dyer-Roeder distance (Dyer & Roeder 1974). Seljak & Holz (1999) 
expanded on this work by considering a mixture of macroscopic and microscopic DM. Unknown to the 
author. Ranch (1991) had previously suggest this use for type la SNe, but had not investigated the shift 
in the mode of the distribution to values below the standard Friedmann-Robertson- Walker luminosity 
distance. Instead he concentrated on the detection of rare high magnification events and concluded that 
many more SNe would be require to detect macroscopic compact objects. In this paper only microscopic 
DM is considered. In this case the DM can be treated as a transparent, massive fluid clumped into 
structures that are much larger than the angular size of the source. 

A few authors have considered how the lensing of type la SNe could be used to probe DM structure. 
Kolatt & Bartelmann (1998) pointed out that the mass sheet degeneracy of galaxy cluster shear maps 
could be broken with SNe. Porciani & Madau (2000) calculated the probability of strong lensing events 
in a CDM universe. They used the Press-Schechter prediction for the mass function of halos and find the 
optical depth for a certain magnification threshold. High magnification events are rare ('^ 2 x 10^"^ chance 
for magnification above 1.3 at z^ = 1) although they could be specifically searched for by monitoring 
high-mass galaxy clusters. The approach discussed in this paper incorporates both strongly and weakly 
lensed SNe to take better advantage of the data and broaden the sensitivity to mass outside halo cores 
hence probing the size scale of the extended DM distribution. The author has previously pointed out that 
lensing by microscopic DM could be detected by an increase in the variance of SN brightnesses at high 
redshift (Metcalf 1999). In addition this increase in dispersion is important for determining the noise in 
cosmological parameter estimates when using SNe at z ~ 1 . The feasibility of measuring lensing via this 
method is potentially hampered by uncertainties in the intrinsic distribution of SN luminosities and its 
possible redshift dependence. In this paper, I seek to eliminate these systematic errors by correlating SN 
brightnesses with foreground galaxies. In this way a null result is expected in the absence of lensing. This 
approach also makes interpretation of the result more straightforward and results in significantly higher 
signal to noise. This work is an extension and improvement an earlier parametric approach summarized in 
Metcalf (1998). 

In section the first order correlation statistic is defined. In section ^ models are constructed for 



estimating and interpreting the contribution of galaxy halos (3.1) and extragalactic clustering (3.2) to the 
lensing signal. Sources of noise are considered in section ^. Some higher order correlation statistics are 
discussed in section ^. The lensing of SNe is compared to the lensing of galaxies in section |6[ The final 
section contains a summary and discussion of future prospects. 



2. Correlating SNe with Foreground Light 

As stated above I seek to reduce as much as possible the influence of SN redshift evolution by 
considering only correlations in foreground light and SN brightness. Using statistics based on these 
correlations also increases the signal to noise and allows for simple interpretations of the results in terms of 
mass-light correlations. Throughout this paper units are used in which the speed of light, c — I. 

To start, consider the weighted flux from foreground galaxies 

J^ ^^w{zi,Zs,yi)ft (1) 

i 

where jji is the proper distance separating the galaxy i from the line-of-sight to the SN, fi is the flux 
received from that galaxy and Zs is the redshift of the source SN. An attempt should be made to correct the 
/i's for extinction and inclination effects as well as k-corrections so that fi accurately reflects the intrinsic 
luminosity of the galaxy. I will choose the weight function to be factorable and normalized as follows 

■w{z,Zs,y) ^wi(z,Zs)w2{y) ,^. 

J^" dz wi{z, Zs) ^l , 27r/g°° dy y W2{y) = 1. 

This gives J-' the units of flux per area. Besides these constraints the weighting function is arbitrary so 
that it can be adjusted to maximize the signal to noise. Later I will make a particular choice of weighting 
function for numerical calculations. 

For simplicity I will consider a flux limited sample of foreground galaxies. In this case the average of J^ 

is 

T= dzwi(z,Zs)—— — ^-r\^) = T~Er Q^ wi(z,Zs) x , . xg / dL L(j){L,z)<i>(L,z) 3 

Jo [^ + z) dz AttHo Jo E{z)dL{zy Jo 

1/9 1 /"^ d?' 

E{z)^[nui + zf + {i-^„,-n^){i + zf + n^\"^ , x(z) = — /-— - (4) 

no Jo -^(.^ J 

where 'q{z) is the number density (number per proper volume) of observed galaxies, dL{z) is the luminosity 
distance, (j)iL,z) is the selection function of the survey and <i>(L, z) is the luminosity function in comoving 
volume. The comoving radial distance is xi^)- The selection function, </>(£, z), will depend on details of the 
particular survey. In the numerical calculations shown later a simple flux limited survey will be assumed 
where (j)(L,z) zero for L < Lf,(z) and one for larger L. In terms of the flux limit Lc{z) = AndLlz)"^ fc- The 
luminosity function of detectable galaxies, (j){L, z)$ (L, z), can be found empirically and independently of 
the SNe. 

A SN is an unresolved near-point source so the only observable effects of gravitational lensing are image 
magnification and the multiplication of images. If DM is microscopic as is assumed here the probability of 
a SN at z 2± 1 having multiple, distinguishable images is very small (< 10"'*). Except for rare cases, SNe 
will be only weakly magnified and demagnified by matter very close to the line-of-sight and no multiple 
images will occur. To second order the magnification, /i = 1 + 2k + 3k^ + 7^. The convergence, k, is defined 
by its relation to the magnification matrix: trA = 2(1 — k) and the shear is 7^ = jA — (1 — k)/|. So in the 
weak lensing limit, ^,7 <?; 1, the magnification is independent of the shear. The convergence, and thus fi in 
this case, depends only on the local mass density (Ricci focusing only) and can be calculated perturbatively 
(Kaiser 1992; Blandford et al. 1991) 

IX f''\i w^ \xf \ w/ ^ 30^go(l + z)g [xjz)] g [xjzs] - x{z)] ,^.. 
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where the density contrast along the line-of-sight is 5{z) — {p — p)/p- The comoving angular size distance 
is g{x) = {Rcurv si'a{x/ Rcurv), X^ RcurvSmh{x/Rcurv)} in the closed, fiat and open cases. The curvature 



scale is R^urv — Ho\J\\ — ^m — ^a|- This can be easily modified in everything that follows to account for 
strong lensing as well as weak. In this notation the luminosity distance is d^lz) = (1 + z)g[x(z)]. 



2.1. Correlations 

Consider the cross-correlation of T with the brightness of the SN, b. At a given redshift the 
magnification is defined with respect to the mean brightness at that redshift. Ideally the average brightness 
of the standard candle, b{z), would be calculated from the SNe themselves. This would make the 
measurement of lensing independent of cosmology and any redshift evolution of the SNe. However this may 
not prove feasible in the short term if there are not enough SNe within a small redshift range to accurately 
determine the average. In this case the cosmology can be fixed so that the redshift dependence of the 
mean brightness is predicted. Determining the cosmological parameters can be effectively done with lower 
redshift SNe which are largely unaffected by lensing and with other observations. The average J-{z) can be 
measured without regard to the presence of background SNe. 



I will consider the correlation function 
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where Abi = bi — b{zi). The sum is over the sample of SNe. If the SNe are used to estimate the average 
b{zi) one must take into account the correlations in the finite sample. The averages of these statistics can 
be expressed in terms of the magnification of the SNe in a redshift range 

{S%;) = {STS^,) . (7) 

where 6T = {T - T) /T. 

To interpret and predict measurements of this statistic, its mean value must be related to the 
distribution of DM and galaxies. For these purposes I define a position-dependent luminosity function by 
writing the probability of a galaxy being within a comoving volume 5V and in a range of luminosity 6L as 

5V6L $(f , L, z) = SV6L $(L, z) [1 + S<P {x, L; %)] . (8) 

The probability of a galaxy being at x is presumably dependent on the density field, S^f, at the point x and 
in surrounding regions. This luminosity function is a function of position and luminosity, but a functional of 
the density field. This functional property expresses the probability that a galaxy's existence and properties 
depend on more than the smoothed density at its position. Nonlocal properties of the surrounding density 
field are probably important. 

In relating these statistics to correlations of galaxies and mass density I will employ an extended 
Limber's approximation which requires that all such correlations are limited to scales much smaller 
than those over which quantities like redshift, the lensing window and the angular size distance change 
appreciably. With this approximation and the luminosity function above, the correlations between 
magnification and J- can be written 



X / dr ^,s,s{\/y^ + r^,L,z). 



where the last radial integral is over proper distance not comoving distance. The luminosity function - 
mass correlation function is defined by 

^*d-(|xi ~X2\,L,z) = {S^{xi,L;Sff)S{x2)) (10) 



where all quantities are evaluated at the epoch corresponding to the redshift z. 



3. Modeling the Signal 

For estimates of the lensing signal, structure in the DM distribution will be modeled as consisting of 
two components. First each galaxy has a DM halo surrounding it whose properties are correlated with 
observable properties of the galaxy. In addition, the galactic halos are clustered into clusters which contain 
additional mass. The lensing contributions from each of these components are considered separately in the 
next two subsections. 



3.1. Galaxy Halos 

Let us make the simplifying assumption that each galaxy sits at the center of a massive halo so that 
(5<I>(x, V) — —1 everywhere but in these locations. In addition, it will be assumed that all the positions and 
properties of the halos are uncorrelated with each other, but directly related to the observable properties 
of the tenant galaxy. These assumptions are probably incorrect in detail, but they will serve well as a 
statistical model whose validity can be tested by the observations. The additional mass contained in groups 
and clusters of galaxies will be handled separately. 

In this galaxy halo model the correlation function is given by 

gj,a(r, L, z) = (5/t(r, L, z) = _ ' ' ' - 1 , (11) 
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where p/i(r, L,z) is the density of the halo which surrounds a galaxy of luminosity L at redshift z and 
T,h{y, L, z) is its surface density. The average is over the orientation and substructure of the halo. 

For the purposes of estimating the signal the surface density of a galactic halo will be modeled by 

v:' 1 .. .. / lV ( l ' 



This profile matches a singular isothermal sphere at small y and then drops away more rapidly beyond r^ 
with an adjustable slope. The central velocity dispersion is related to the luminosity by the TuUy-Fisher 
relation, /? ~ 1/4. The scale size, Tc, is related to the luminosity by a less certain relation. In what follows 
a = 1/2 which in the case oi P — 1/4 gives a constant mass to light ratio for the isothermal part of 
the halos. In addition this makes Tc proportional to the radius of the visible galaxy in an idealized case 
where the surface brightness is a constant. The luminosity function is taken to be a Schechter function 
^{xL^,) = <I>*.T^^°^e^^ with L* = 10^'^ H^'^Lq and independent of redshift. The weighting function will be 

W,iz,Zs)-W,iz,Z,)dUzr , W,iy)=l^^^^ Rrmn<y<R ^^^^ 

y y > R , y < Rrmn 

This is not necessarily the optimal weighting, but it does seem to give high signal to noise without 
prejudicing the result by assuming a specific halo profile. With this weighting function S^ (R) is directly 
related to the average surface density in the annulus Rmin < y < R surrounding a galaxy. 

In the above the weak lensing approximation has been used. More generally SjrJR) is related 
to the average magnification within the annulus. However the probability of forming an image is 



not uniformly distributed on the image plane due to magnification bias. The result of this is that 
S^^{R) oc {fJ.)sn — 1 = (m~^)x ^ ~ 1 where (. . .)x is the average over the annulus Rmin < y < R- Weak 
lensing will be assumed throughout this paper so this distinction in the interpretation of S'^ (R) will not 
be important. 



3.2. Extragalactic Clustering 

It is well established that in the central regions of galaxy clusters most of the mass is not associated 
with individual galaxies but resides in intergalactic space. This is confirmed by observations of gravitational 
lensing using galaxy shear (Natarajan et al. 1998), the Sunyaev-Zel'dovich effect and X-ray emission. This 
extra cluster mass can have significant lensing effects and conversely any information that can be gained 
through lensing on the structure of clusters would contribute to our understanding of cluster formation and 
evolution. 

To model the effects of extragalactic mass - mass outside of any galactic halo - I treat the galaxy halos 
as being embedded in larger halos with a universal profile but of varying size and mass. Only the component 
of extragalactic mass that resides in virialized halos will be considered because any other component will 
not contribute appreciably to the lensing considered here. Because 5/^ oc S to first order the contribution of 
the extragalactic mass to lensing can then be calculated separately and simply added to the contribution 
from galactic halos. This will not be true for strongly lensed SNe, but these are comparatively rare and 
will not be considered here. I will call the extragalactic contribution the cluster contribution although 
it will include extragalactic halos with masses well below what is normally called a galaxy cluster. The 
division between galactic and cluster mass is somewhat arbitrary in any particular case, but it does have 
a meaning in a statistical sense. By comparing cluster galaxies with the 'average' galaxy the extra mass 
can be attributed to the cluster or group. This separation is more direct and less model dependent than 
it is for shear measurements (see Natarajan, Kneib, & Small 1999; Natarajan et al. 1998) owing to the 
magnification being directly proportional to the surface density. 

The relative mass density of an extragalactic halo is written 5c{x^ M) and its projection onto two 
dimensions is SYic{y, M). If the density of galaxies is assumed to be proportional to the overdensity the 
correlation on the sky between galaxies and surface density in halos with masses between M and M + dM 
is given by 



{y,M) - f^{M)fg{M) f^fd^x 5Y.,{x,M)5Y.,{x-y,M)l j d^'x S,{x,M) (15) 

«i'fi«^^*^A(MI*E.(M/)P (16) 
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where fg{M) is the fraction of galaxies in halos within this mass range and fY.{M) is the fraction of the 
cluster surface density that is not in any galaxy halo. The first integral in (|l5| ) is an average over the 
direction of y or cquivalcntly the orientation of the halo. It is often more convenient to calculate the 
correlation function from the Fourier transform of the halo profile, 5Y,c{k^ M). If the weight function W2{y) 
is given by (nj) with Rmin = 0, a top hat, the relevant quantity is the average correlation function within a 
circle of radius R, 

^(^^'^) = ^p^wn^fn /* -^i^^^) m{k,M)\\ (17) 

7ri?dEc(0, M) J 

The J's are Bessel functions. 

Computationally the problem can be simplified when the halo profile obeys the scaling law 

<5Sc(fc, M) = 5,{M)r,i{Mf5i:, {rci{M)k) ; i%{y, M) = f^{M)fg{M)5o{M)r,i{M) i% [jjj^yj (18) 



where rci{M) is the scale length of the halo and 5o{M) is a density scale. Many popular profiles are of this 
generic type. The scaling implies a contribution to S^ from extragalactic halos of 

{S%)^^Y.T- r dzw,{z,z,)^l±^^^^ r dLL^{L,zW\L,z) (19) 

/•CO / n 

X / dM MM, z)/,(M, z)6oiM, z)r,i{M, z) i% 

This scaling allows one to calculate G^{x) once and then integrate over halo mass and redshift. The mass 
Mmin is the mass of the smallest halo that contains multiple galaxies. The luminosity function is given a 
superscript because there is a known segregation of galaxy types between clusters and the field. However 
in the calculations to follow I take the cluster galaxies to have the same luminosity function and halo 
distribution as the galaxies in the field. 

For signal estimates I will assume the extragalactic clusters are of the Navarro, Frenk and White 
(NFW) form (Navarro, Frenk, & White 1997, 1996): 

^ciir)- , , .,f° , ., . (20) 

{r/rcij (1 + r/rrj) 

The average galaxy-surface density correlation function, ^^^ (x) , in this model is plotted in figure || This 
function is effectively constant for x ;$, 0.1 and zero for x ^ 10. 

To calculate dl3) some assumptions concerning the population of halos and their galaxy content must 
be made. I will consider fiat CDM models with scale free {n = 1) initial perturbations. The parameters 
So{M, z) and rc{M, z) are related by re = r2oo{M)/c and So — gjFTTyC^/ [ln(l + c) — c/(l + c)] where 
''200 (M) is the radius within which the mean density is 200 times the background. To simplify the procedure 
of calculating these quantities the simple formula c ~ c*(M/M*(z))~°-^ is used. The nonlinear mass scale 
is defined by the condition Ao(M*,z) = 5c{z) where the critical density is 6c — 1.6865 rJm(z)°°°^^ for flat 
models. These scalings reasonably reproduce the results of Navarro et al. (1996). The normalizations for 
the ACDM model considered here (fi^ = 0.3, Oa = 0.7) are CTg = 1-14 and c* = 7.9. For standard CDM 
{flm = 1), CTg = 0.6 and c* = 10. The fluctuations on 8h~^ Mpc scales, as, is consistent with observed 
galaxy abundances in both models (Viana & Liddle 1996). 

In populating the clusters with galaxies it is assumed that the fraction of galaxies in clusters within 
a given mass range is equal to the fraction of total mass which is contained in those clusters. The mass 
fraction can be calculated using the Press-Schechter prediction for the mass function of halos. This results 
in 

The galaxy halos are statistically the same in and out of clusters. The fraction of the mass in the cluster 
but not in a galactic halo, /s(M, z), will not depend on M in this case. The assumption that the galaxy 
content of a cluster is proportional to mass results in f^{z) = 1 — pg{z)/'p[z) where Pg{z) is the total 
mass content of galaxies and their halos. With the static halo model considered here /s is independent of 
redshift. Typically /s ~ 0.4. This number would increase if cluster galaxies are tidally truncated as they 
undoubtedly are in reality. 

Figure shows the expected size of S^ [R) for SNe at z = 1 where R is the radius of the weighting 
function in w%. On top are curves representing the contribution of galaxy halos. In each case an artiflcial 
inner cutoff, Rmin , is made so that SNe whose line-of-sight comes within this distance of a galaxy are 
excluded in order to avoid problems with obscuration. This is an adjustable parameter that is applied 
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Fig. 1. — The galaxy-surface density correlation function averaged over a disk, Coe(^)- ^^'^ cxtragalactic 
halo is taken to be of the NFW form (equation EOl) . 
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Fig. 2. — The correlation of light and magnification SjrJR) for SNe a,t z = 1. The contributions from galaxy 
halos and extragalactic halos are plotted separately. The total <S|r„ is the sum of the two components. The 
galactic halo curves are the steeper ones starting in the upper left corner. Two different models which 
differ in the logarithmic slope of the surface density at large galactic radii are shown. For each model two 
inner cutoffs are considered ~ Rmin — 10, 20 kpc. All SNe within Rmin of any galaxy are excluded. The 
two curves that flatten out on the left are the contributions from clusters in the CDM model. The one 
marked fg{z) = 1 represents the case where only SNe behind M > 10^"' M0 clusters are selected. The 
crisscrossed region represents the uncorrelated noise contributed by uncertainties in SN peak luminosities 
{csn = 0.16 mag) and shot noise in the galaxy counts for the case of 250 SNe (see section ^). The signal 
to noise ratio at i? = 200 kpc is ~ 4-y/iVs„/250. The limiting magnitude of foreground galaxies is m = 23. 
Additional parameters are printed on the plot. 



when analyzing the data. SjrJR) is then proportional to the average surface density between Rmin and 
R. Figure ^ shows that S'^ (R) is relatively independent of Rmin when R ^ Rmin- The inner part of 
these curves, beyond the Rmim depend only on the isothermal part of the halo model and thus only on 
the TuUy-Fisher relation. This could in fact be a way of determining if a TuUy-Fisher-Faber-Jackson type 
relation applies to halos. Two different halo models are shown with 7 = 2 and 1.1 (halo model is given 
in (13)). These models start to diverge from each other at large R (r* = 500 kpc in these models), but 
differentiating between them is likely to require a large number of SNe both because of the small difference 
and the contribution of cluster mass which starts to become significant around the same radii. It is clear in 
both models that SjrJR) deviates from isothermal as R gets larger. A measurement of the halo size scale, 
r*, should be possible. 

The cluster contribution is represented in figure by the lower curves which flattens out at small 
R. The lower of these curves is the contribution expected for a random sample of SNe. In the present 
model this contribution is fairly small relative to the noise. This is largely because of the small number 
of galaxies that reside in large clusters. For these calculations the minimum mass of a cluster is taken to 
be Mmin = 3 X 10^^ M0, but it is the larger clusters which give most of the signal. All else being fixed, 
the extragalactic halo component will be proportional to fg (M) - the fractional galaxy content of clusters. 
In this way it is a measure of bias. The noise for a sample of 250 SNe is represented in figure g by the 
crisscrossed region. With this number of SNe and the level of noise assumed here the signal is ~ 4.5 x cr. a 
detection at the 4.5 cr at i? ~ 0.1 Mpc. The noise is discussed further in section^. 

The relative size of the cluster contribution can be increased by selecting SNe - or searching for SNe - 
which occur behind identified clusters. The upper cluster curve in figure marked /s = 1 shows the level 
of signal if all the SN sight lines are within ~ 10 Mpc of the center of M > 10^"* M© clusters. Now the 
contribution from extragalactic material is comparable or greater than the contribution from galaxy halos at 
R ^ 1 Mpc. In actuality halos within a cluster should be tidally truncated which will make the halo curve 
fall more rapidly with R and the cluster curve will he higher due to the redistributed matter. It would be 
interesting to study the tidal striping process in this way. It is strongly dependent on the formation history 
of clusters. Of course the cluster mass threshold could be increased so that the cluster lensing contribution 
is comparatively larger. In this way the substructure in galaxy clusters can be measured very directly. 

Figure H shows the same SjrJR) plot with the background cosmology changed to Qm = 1- The galaxy 
halos arc the same as in figure 0. This galaxy halo signal is almost independent of cosmology. Perhaps 
unexpectedly A does not have a large effect even when the galaxy halo model is held fixed. This is a result of 
using correlations between light and magnification rather than a statistic based purely on the magnification 
probability distribution such as the optical depth. The magnification probability is a stronger function of 
the angular size distance and the density of halos while tS|r„ is normalized to fit the observed density of 
galaxies. The cluster component is a little bit bigger in the high Qm models when the power spectrum is 
normalized to observed cluster abundances. 

The modeling of extragalactic halos done here is somewhat crude. There are many ways in which it 
could fail to represent reality in detail. The density of galaxies within a halo is not necessarily proportional 
to the halo density and the total number of galaxies in a halo may not be proportional to the halo mass. It 
is already known that galaxies of different morphologies, colors and luminosities are biased with respect to 
each other. This would seem to indicate that these assumptions are not entirely correct. In addition, large 
clusters tend to have a central elliptical galaxy which has not been included here. This would increase the 
statistical weight of the cluster contribution. More sophisticated modeling is certainly possible when the 
data warrants it. 
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4. Noise in {S'jr ) 

There are several different sources of noise tlrat will affect the precision and accuracy of any {S"^ [R)) 
measurement. First there is noise coming from the dispersion in SN luminosities and the random fluctuations 
in foreground light. This noise is uncorrelated with the lensing signal. The lensing structures themselves 
will also contribute to the noise through sample variance. The density within radius R is concentrated 
toward the center of the halo. The more concentrated it is the less good an estimate of the average from 
a finite sample is likely to be. This problem is compounded if the DM halos are not smooth but clumped 
into substructures. In addition errors in the redshifts of foreground galaxies will add noise. There are also 
observational errors and biases arising from the light curve correction to the SN peak luminosity, the host 
galaxy subtraction and the detection efficiency of the SN search. For lack of detailed information on future 
observations I will not attempt to address the last categories of noise, but I will estimate the importance of 
the other sources of error. 



The variance in the correlation can be written 

r . . . /A7;\^ . 1 

(22) 



(['^JPmJ ) ^ {^Fm) — =2 

NsnT 






where A6(nm,riA) is the difference between the true average brightness and the one derived using the 
parameters; Vtm,^A, etc. The variance in SN brightnesses caused by intrinsic differences and observational 
errors is tr^. It is assumed that the error in T is small enough that it is unimportant. The last term 
represents the sample variance produced by the lensing itself. All of these quantities will be redshift 
dependent so in a real sample they will be different for each SN. 

Errors in h{^rm ^a) should always be a subdominant. This is because the cosmology will be constrained 
by all of the SNe observed, including those at lower redshift, so it is expected that A6 '^ <y^- Of course 
there is always the possibility that one has chosen the wrong parameters to describe the real universe. 
With enough SNe one could directly calculate hiyim^^h) within redshift bins. This would make the result 
independent of cosmology at the expense of increasing the h{^„n ^a) term. 

The variance of JT is best estimated directly from data. However, to assess the noise levels expected in 
a hypothetical experiment it is useful to make a semi-empirical model for the variance a'jr. If it is assumed 
that the internal properties of different galaxies are uncorrelated the variance is 

2^ r^ ^^ wi{z,z,Y{l + zY 
Ho J a ^ [AT:dL{z)^f E{z 

+ ^flo^dLLcl^{L,zML,z)j;;^dL'U<j>{L',zML\z)fd^y,fd^y,w,{y,)w,{y2) 

(24) 



d^ : ::\.\2l:\ \ I dLl^d^mzM^z)! dyyw^iyf (23) 



jr^dr^,,[^/A^^T^,z)} 



The first term (|23| ) is from shot noise in galaxy number counts. This expression can be derived by 
considering the variance in the galaxy number count within a region of space small enough to contain at 
most one galaxy. The probability of a galaxy being in this volume is given by equation m) . Then a Poisson 
distribution implies 

(<^{L)^{L')S^{x, L)6^{x', L')) ^ <^{L)6{L - L')S^{x - f')/(l + zf (25) 



where the positions x are measured in proper coordinates. Term ( |23D follows from this and is plotted in 
figure ^ for the low J^a flat model. 

The second term in ctJt, (p4[), is from the clustering of galaxies. It is known from observations that 
the correlation function of galaxies is well approximated by £,gg{x) = {x/ro)~^ (7 = 1.77 ± 0.004 and 
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To = 5.4 ± Ih ^ Mpc) and independent of redshift if x is the comoving distance. With this level of clustering 
(7^ is dominated by shot noise for relevant range of R so clustering contribution will be ignored. 

In addition to shot noise the estimate of J- is limited by errors in the estimated redshifts of the 
foreground galaxies. These redshifts need to be estimated either photometrically or spectroscopically. 
This can cause errors in two ways. If one is using the weighting function (04) errors in redshift can cause 
'Wi{zi, Zs) to be misestimated and it can cause a galaxy that is outside the volume allowed by W2{y) to be 
misidentified as inside that volume and vice versa. If necessary the latter problem can be avoided entirely 
by using an angular cutoff in W2{y) instead of a proper distance cutoff as in (|lj). This is done at the 
expense of being able to interpret {S'^ (R)) as a measure of the average surface density within distance R 
of a galaxy. An estimate of the uncertainty in J^ caused by errors in wi(z, z^) can be found by replacing 
wi{z, Zs) in ( [2^ ) with 7]4pjy(l + ^)^^{^7 ^s) where Sz/{1 + z) is considered constant. This is also plotted in 
figure 0. In the case considered here (JredsMft — 10J2;/(1 + z)ashotnoise- At this time photometric redshifts 
have an accuracy of 6z ~ 0.1 — 0.2 for the relevant redshift range (Bolzonella, Miralles, & Pello' 2000). If 
photometric redshifts are used this noise is comparable to the shot noise contribution. With spectroscopic 
redshifts this noise would be unimportant. In figures and only the shot noise times the SN luminosity 
noise is shown. 

The size of the sample variance depends on the structure of the halos. If the halos are very centrally 
concentrated, asymmetric and/or clumpy the average S^ will strongly depend on a small number of SNe 
which happen to be observed through high density region. In this case S^ may converge to its mean value 
very slowly and a prohibitively large number of SNe may be required to estimate it. Ignoring correlations 
between halos the sample variance can be written 

jr^^ Jo [47rdi(z)2]>(z)»-2 

X J dL L^(b{L, z)$(L, z) J d?y w^ [yfB'^'-^iL, z, Zs, y) 
with n — A (the n 7^ 4 will be useful later). 

In the weak Icnsing limit B"{L, z, Zg, y) — (S(i, J/)")- If the halo structure is perfectly correlated with 
the galaxy luminosity B^^{L, z, Zs,y) = Y,{L^yY^. There is undoubtedly some scatter in the luminosity- 
surface density relation which will increase Vjr^. The scatter is not known, but if we take the scatter in 
the TuUy-Fisher relation as a guide Vj^^ will be about 10% larger than the perfect correlation value. This 
should be added to the lower solid curve in figure 0. The sample variance would then be smaller than the 
first term in (^). However the scatter could well be larger than it is for the TuUy-Fisher relation in which 
case one might hope to measure it with higher order statistics. This is discussed in section p.l| . 

Another possibility is that the there is a significant amount of substructure in the dark matter halos. 
In this case I will resort to an approximation for _B"(z, Zg, y) to simplify calculations. In order to preserve 
the average angular size distance the average magnification at a distance y must be the same with or 
without substructure. A model for substructure lensing should respect this constraint and take into 
account the strong lensing that may be important if subclumps have singular cores. I consider the case 
where all the mass in the halo is contained in randomly placed clumps each with average surface density 
Tisub- The probability of a clump being in the line-of-sight is then p = Y,{y)/T,sub- Outside the clumps 
the magnification is zero (it depends on shear only at second order). When the halo density is large the 
possibility of two clumps eclipsing the line-of-sight must be taken into account. In the weak lensing limit 
B"-{z, Zs,y) ~ p{l + 2p) (Ssub) where only the terms of highest order in T,sub have been retained. To take 
strong lensing into account I renormalize this by ((5/i")/((5/i)". The result is 

{E(2/)" , no substructure 

|gfPsr;S(L, ,) (1 + 2E- E(L, ,)) , substructure ^''^ 
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If a magnification cutoff is imposed to make ((5/i") converge, ((5/i)" is still calculated without this cutoff so 
that the proper normalization is retained. 

To model the magnification within a subclump I use the solution for an isothermal sphere cutoff 
at an outer radius, Rsub- The characteristic length scale for such a lens is Xo = 4:TrvlDiDis/Ds. The 
magnification, including both images, is then a function of only the rescaled impact parameter x = y/xo'- 





- (2-x)/a; 


, X <1 


Sfi = < 


l/x 


, 1 < a; < a; 







, X > Xm 



(28) 

The second moment is logarithmically divergent so a maximum magnification cutoff must be imposed to 
calculate (p9). This cutoff can be set high enough that the probability of there being a SN in the observed 
sample with a magnification this high is very small. The divergence of the moments is of course not real. 
In actuality the finite size of the source puts a limit on the magnification. 

The size of dark subclumps is not very well constrained by observations, but they do appear in 
simulations. Moore et al. (1999a) find a large number of subclumps in their CDM simulations - many more 
then are accounted for by visible dwarf galaxies. For these simple estimates I take the velocity dispersion 
to be Vc = 15 km s^^ and the size to be 100 pc. This is an extreme model since all the mass is in rather 
large subclumps, but it has the advantage of simplicity and gives a conservative estimate for the noise. The 
sample variance, VJr^{R), with and without this substructure is plotted in figure |j. The cutoff is set at 
Sji > 10. In this model there is about a 4 x 10"'* chance of a SN at z = 1 being magnified above this limit. 

With the variance expected in the corrected SN luminosities the uncorrelated noise, the first term 



in (22), generally dominates over the sample variance. At small R the sample variance may be comparable 
to the uncorrelated noise, but to some extent this is a result of the model used here. One would expect the 
subclumps to break up when the density is high which is not reflected in the model. Another reason to 
think that this is an overestimate of the sample variance is that the density profiles of real CDM clumps 
are probably not as steep as p ex r~^ at small r. It has recently been claimed that p oc r~^'^ is a better 
approximation (Moore et al. 1999b) and the NFW model has p oc r^^. In these cases T'^ would be smaller. 

There is one other outstanding problem that could still cause trouble for this measurement. It has 
been suggested that an exceptional kind of dust could account for the observed faintness of high redshift 
SNe. The dust grains would need to be large enough to avoid the constraints on reddening (Aguirre 1999; 
Perlmutter et al. 1999) and quite uniformly distributed in the intergalactic medium (Croft et al. 2000). 
This is not consistent with the properties of known dust, but it does remain a possibility. The distribution 
of this dust would presumably be correlated with galaxies so it would show up as an anti-correlation 
between foreground light and SN brightness. Looking for correlations might actually be the best way of 
testing this hypothesis as well. 

In summary it seems promising that (>Sjr„) could be measured with ~ 250 at z ~ 1. If the SNe are 
observed with a future satellite it is likely that the noise will be smaller than is assumed here because in 
this case agn will be smaller at high redshift. 



5. Higher Order Correlations 

The above discussion of sample variance leads to the question of whether dark subclumps in galaxy 
halos might be detectable through the gravitational lensing of SNe. Such a detection would have important 
implications for structure formation theories. Recent high resolution N-body simulations have revealed that 
the CDM models have several apparent problems when compared to observations. One such problem is 
that galaxy halos in simulations have much more substructure or smaller clumps within them than observed 
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Fig. 3. — The same as in figure except flm — 1, h — 0.5 and as = 0.6. 
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Fig. 4. — The major sources of noise for the statistics S^ and 7^^„. These are all calculated with the 
parameters used for figure ||. The dotted curve is the noise just from errors in redshift times (1 + z)^/(5z^ 
which is regarded as a constant. The dot-dashed curve is the noise from shot noise in the galaxy counts. 
For each pair of P^^ curves the lower one is without substructure and the higher one has substructure. The 
substructure clumps, as described in the text, are truncated isothermal spheres with radius Rsut = 0.1 kpc 
and velocity dispersion Vc = 15 km s^^. The magnification cutoff for P^ and 7-"^ is set at Sfi — 10. All 
quantities are dimensionless. 
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galaxies have observable dwarf satellites (Moore et al. 1999a). These subclumps form early and fall into the 
larger halos without being destroyed by tidal stripping and collisions. Subclumps of dwarf galaxy size seem 
to be over produced by about a factor of 50 in halos that should host Milky Way-like galaxies. Considering 
the successes of the CDM model it seems reasonable to entertain the possibility that these subclumps exist, 
but are dark either because they never formed stars or because their star formation was shutoff early.R 

In this section I investigate the two third-order correlation functions of 6fi and JF and find that one is 
sensitive to this substructure. The other contains additional information on the clustering of dark matter 
on extragalactic scales. 



5.1. Light-magnification-magnification Correlation 



In direct analogy to S'^ I define the statistic 7?.jr 



^i..-^E^V ' (T^h) = i^^S^') (29) 



3 ■' J "j 



Like 5|r„ the average of this statistic will be zero in the absence of lensing. 
Assuming weak lensing the average value of 7?.^^ is 

\4 7P/'.,\TT/ /.v\2 



(^^p)-^ = :tt^ >^ / dzwi{z,Zs) T~ry2 / dL L(f>{L,z)^{L,z) dy yw2{y) (30) 

^JVsn ^ Jg O-LlZ) Jo Jo 

dridr2 C$55 f 2/, yy'^ + rl, \J y'^ 



OC /"OC 



-^r^,L,z 



oo >/ — oo 



C*55(xi,a;2,a;3,L,z) = {5<^{xi, L; 8^)8{x2)5{x'i)) 

'^%ii '^^ related to how fluctuations in mass density, instead of just the mass density, are correlated with 
light. Taking each galaxy halo to be independent and thin results in the third order analog of equation (O) 



OC /"OO 



oo -/ — oo 



dr^dr2 C^ss \v, \l v' + rl\l V^ + rlL^z] = ^^j;!^;^' ^jP^ ■ (31) 



The second moment of the surface density will differ from simply the square of the average halo profile 
because of scatter in the L-E relation, asymmetries and substructure. Asymmetries may increase the 
moment by a small factor. Scatter L-E relation could potentially be larger - perhaps a factor of a few. 
Substructure could potentially have a largest effect. I use the substructure model discussed in section to 
estimate how sensitive 7?.^^ is to substructure. The resulting calculation is plotted in figure ||. 

The variance in 7?.jr„ can be written 

( [7^^„] ') - {n%j^ = ^ {cy^^^,l + s%^ q%^ + v%^l + v%y, + v%^) (32) 

where P^ is defined in (Effl and 



Q%^^-4=\ dzw,{z,z, y \ rJ m{z,zs)^ I dLL<l>{L,zML) j d^y W2{y)B\L,z,Zs,y). (33) 



—= dzwi{z,zs) -r-TT 

inT Jo dL{z) 

'^It is interesting to note that if this is the case the present microlensing results do not put very strong constrains on the 
fraction of the halo in MACHOs (see Metcalf & Silk 1996). 
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Fig. 5. — The three sohd hnes are the second order correlation function 7?.^„. The lowest of these is the 
result with no halo substructure and perfect correlation of the halo with the luminosity of its resident galaxy. 
The dotted curve is the noise in this case. The scatter would need to increase 7?.^^ by a factor of several 
for it to be measurable. The other curves include substructure in the halos. The top solid curve is the result 
with 6)1 > 10 excluded. All the other curves are for 6)1 > 4 excluded. The inner crossed region represents 
the noise without sample variance. The dashed lines are the one-c noise with sample variances. These are 
all calculated for 800 SNc at z = 1. 
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The approximated value of B'^{L, z, Zs,y) is given in (pTf). The moments of the distribution of intrinsic 
SN luminosities are /i" This distribution will be approximated as Gaussian in which case /i^ = and 
n't — 3(7^. In reality each of these factors will be different for each SN because of differences in redshifts 
and observations so an average over SNe will be done. 

Unlike in the case of S^ the noise in TZ^^ is usually dominated by sample variance rather than SN 
luminosity errors and galaxy count noise. As a result both the signal and noise are strongly dependent on 
the structure of the subclumps. This is a result of the strongly peaked central density of the subclumps. 
Specifically the Vjr,^ term in (|3^) tends to dominate. If the cutoff in 6^1 is high the sample variance will 
swamp the signal. As can be seen in figure a if there is no substructure (7^jr„} is too small to be measured 
with 800 SNe at z ~ 1 because of sample variance. Substructure increases the signal, but also the sample 
variance. A magnification cutoff can be used to reduce the sample variance without reducing the signal as 
significantly. The cutoff is a censoring of the data in order to improve the convergence of the remaining 
subsample. The high magnification events are rare so if the cutoff is at high magnification the subsample is 
likely to contain all the data. The subsample can be viewed as sampling the area on the surface of the halo 
where the magnification is below S^max- By adjusting Sfimax we can choose to have a more precise result 
relating to a restricted area or a less precise result relating to more of the area. 

The sample variance may not be as big a problem in reality as it is in this model of substructure. 
Pure CDM simulations indicate that the density in the clump cores is not as steep a function of radius as 
in the singular isothermal sphere model (Navarro et al. 1996; Moore et al. 1999b). The density is more 
likely to go as r~^-^ or r~^ rather than r~^. This reduces (7^|r„), but it reduces the sample variance more 
significantly. On the other side, galaxy halos are not likely to be made entirely of such large subclumps. 
Moore et al. (1999a) find that subclumps in their simulations have a mass distribution dn/dm oc m^'^ down 
to their resolution limit. Smaller or fewer subclumps would reduce the signal by approximately the covering 
fraction of the clumps. 



5.2. Light-light-magnification correlation 

Finally, for completeness, there is the light-light-magnification correlation function 



t3 = 



J-'fl 



^^A^ - {lh)-{^^'^^^) (34) 






i^h)^' = it^^Jo'd^ u;i(z,z,)^ y^+-);^'i;; /„°°dL Lcp{L,zML,z)j^dL' L' cl^{L' , z'^L' , z) 



X !-^d^f ST^d'^r' W2{ri_)w2{r'^_) C$<i.5 (M, \r'\ 



\r — r 



.L,L\z 



(35) 



C,^<s,s{xi,X2,xz,L,L' ,z)^ {5^{xi,L;S^)S^{x2,L';S^)S{x3)) (36) 

The integrals over the third order correlation function Ci^s{xi, X2,X3, L, L' , z) can be expressed as the sum 
of two terms. One is a shot noise term resulting from galaxy halos taken individually. Like {Sjr^) this term 
is proportional to the galaxy-surface-density correlation function. Then there is a term from correlations 
between different galaxies which is proportional to the galaxy-galaxy-surface-density correlation function. 
This term falls off less rapidly with R so that at large radii {I'jr^ is related to the average surface density 
associated with multiple galaxies. 



It was shown in section 3.2 that by selecting (or searching for) SNe that are seen through galaxy 
clusters SjrJR) can be used as a probe of cluster matter. The selection process would likely have a 
rather high cluster mass threshold. To probe smaller mass clusters and halos a more systematic approach 
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might be useful. The statistic IjrJR) provides a more direct, although somewhat noisier, way of probing 
extragalactic halo structure. This will not be feasible until a considerable number SNe have been observed 
so no detailed calculations of (X|r ) will be attempted here. 



6. Comparison with the lensing of galaxies and quasars 

The gravitational lensing of galaxies is now an established field of observational astronomy. It is 
appropriate that the lensing of SNe be compared with that of galaxies to argue for their independence 
and complementarity. Weak lensing can be detected either through the distortion of galaxy images or 
variations in their number counts. The former is a measure of the shear field and the latter is a measure of 
the magnification. Galaxy-galaxy lensing, the technique most similar to the one discussed here, is where 
the image distortion of background galaxies is correlated with their projected distances from foreground 
galaxies. This results in a measurement of the average shear as a function of distance from a galaxy. Such 
measurements have been done or attempted by several groups (Fischer et al. 1999; Hudson et al. 1998; 
Brainerd et al. 1996; Griffiths et al. 1996; Tyson et al. 1984). The tangential shear averaged around 
any circle on the sky of radius 9 is given by the remarkably simply relation Jt{0) = k(< 9) — 'K{9) where 
k(9) = T,{9)/Tjc{z,Zs) is the convergence. The average convergence within the circle is k(< 9) and the 
average on the circle is k(0). In the weak lensing limit Sfi{9) — 2k{9) so for a singular isothermal sphere 
Jt{9) = 'k{9) and 6ii{9) ~ 2jf{9). If the mass profile falls off more rapidly than isothermal, as it must for 
the mass to be finite, 6^{9) will fall off' more steeply than 7((^). For this reason measuring i5/i(r) with SNe 
might be a better method for finding the size of galaxy halos. Galaxy-galaxy lensing has not been able to 
put any clear constraint on the size of halos. The transition between galaxy halos and intercluster material 
might also be more distinct with SNe. 

The shear and magnification measurements could also be combined; galaxy-galaxy lensing constraining 
mostly the mass scale of halos and SN lensing the size scale. In this way degeneracies and systematic errors 
could be significantly reduced. The signal to noise per source is actually better for SNe than for galaxies. 
The rms ellipticity of galaxies is (e^)^/^ ~ 0.3 so for isothermal lenses roughly 10-25 times the number of 
background galaxies as SNe are needed to attain the same signal to noise. This factor brings the number of 
require SNe into the range that is already being considered for measuring cosmological parameters. 



As discussed in section 5.1 there is the possibility of measuring substructure in galactic halos with SNe. 
Because the image sizes of the background galaxies are generally much larger than the size of subclumps 
galaxy-galaxy lensing will not be able to measure any realistic level of substructure. In addition, the 
contributions to the shear from the halo as a whole and the subclumps are incoherent. In contrast the 
lensing of SNe is sensitive to substructure masses all the way down to a small fraction of a solar mass. 

Another use for the lensing of galaxies is to actually map the shear field. Because each galaxy has an 
intrinsic ellipticity the shear can only be measured in cells which contain multiple galaxies. This limits the 
angular resolution on which the shear field can be mapped. The same is true for magnification bias where 
the number density of background galaxies is related to the magnification. The lensing of SNe is unlikely 
to ever be useful for this kind of study because the number density on the sky will always be significantly 
lower than the number density of galaxies. However, SNe may be used to remove the mass sheet degeneracy 
that exists in shear maps of galaxy clusters (Kolatt & Bartelmann 1998). 

Quasars have been used as sources in lensing studies for years. They are similar to SNe in that they 
are nearly point sources, but they are not standard candles. As a result lensing is only detectable when 
multiple images are formed. This is a rare occurrence and tells use more about the central regions of halos 
than their outskirts. Monitoring the light curves of multiple images does seem to be a promising method 
for detecting compact objects. Subclumps the size of dwarf galaxies however will have time scales that are 
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to long to be detected in this way. 



7. Discussion 

It has been shown that gravitational lensing can be measured by correlating SN brightnesses with the 
density of foreground galaxies. The result is a measure of the correlation of mass with light, the bias, on 
^ Mpc scales. The size and mass scales of galaxy halos should be measurable with a few hundred SNe 
at z ~ 1 . The structure of galaxy clusters can also be probed by selecting SNe that are viewed through 
clusters. Higher order correlations can be used to probe substructure in galaxy halos. In this case the 
convergence of the result is strongly dependent on the form of the substructure that is present. 

The prospects for making these measurements seem very good. The numbers of SNe and the noise 
levels discussed in this paper are attainable, even conservative. The proposed VISTA telescope is expected 
to be able to find hundreds of SNe at redshifts ;^ l.|j The proposed SNAPSAT satellite could find thousands 
up to z ~ l-7.n At high redshift the major source of uncertainty in a SN luminosity comes from subtracting 
the light of the host galaxy. The superior resolution of a satellite could reduce the variance in the corrected 
peak luminosities of high redshift SNe to that of local SNe, Am ~ 0.12 mag or better. This would be a 
significant improvement on the Am ~ 0.16 mag assumed throughout this paper. Gravitational lensing adds 
strongly to the motivations for these projects. 

I would like to thank D. Mortlock and P. Natarajan for very helpful comments on this paper. 
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